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Surface transport of inertial particles is investigated by means of the perturbative ap- 
proach, introduced by Maxey (J. Fluid Mech. 174, 441 (1987)), which is valid in the case 
the deflections induced on the particle trajectories by the fluid flow can be considered 
small. We consider a class of compressible random velocity fields, in which the effect of 
recirculations is modelled by an oscillatory component in the Eulerian time correlation 
profile. The main issue we address here is whether fluid velocity fluctuations, in par- 
ticular the effect of recirculation, may produce nontrivial corrections to the streaming 
particle velocity. Our result is that a small (large) degree of recirculation is associated 
with a decrease (increase) of streaming with respect to a quiescent fluid. The presence of 
this effect is confirmed numerically, away from the perturbative limit. Our approach also 
allows us to calculate the explicit expression for the eddy diffusivity, and to compare the 
efficiency of diffusive and ballistic transport. 



1. Introduction 

Particle transport in laminar/turbulent flows is a problem of major importance in a 
variety of domains ranging from astrophysics to geophysics. For neutral (i.e. having the 
same density of the surrounding fluid) particles in incompressible flows the main quantity 
of interest is typically the rate at which the flow transports the scalar, e.g. a pollutant. 
For times large compared to those characteristic of the flow field, transport is diffusive 
and is characterized by effective diffusivities (the so called eddy diffusivities) which in- 
corporate all the non trivial effects played by the small- scale velocity on the asymptotic 
large-scale transport ( Biferale et a/.l 19951 Mazzinolll997h . Eddy diffusivities are expected 
on the b asis of central limit arg uments and can be calculated by means of asymptotic 
methods ( Bender fc Orszaj|l978 ). There are however situations where some, or all, of the 
hypotheses of the central-limit theorem break down with the result that in the asymptotic 
limit particles do not perform a Brownian m otion and anomalous diffusion is observed 
( Castiglione et al. 19991 : Andersen et al. 2Qq3)- 

For compressible flows, large-scale transport is still controlled by eddy diffusivities but 
now transport rates are enhanced or depleted depending on the detailed structure of the 
velocity field (|Vergassola fc Avellanedalll997l ). 

For particles much heavier than the surrounding fluid, large-scal e transport is st i ll con- 
troUed by eddy diffusivities, which have been calculated in iPavliotis fc Stuart ( 2005h 



2 



M. Martins Afonso, A. Mazzino and P. Olla 



exploiting asymptotic methods in the limit of small inertia. 

Unlike what happens for neutral particles, an external force field can dramatically change 
transport properties. This is for instance the case when gravity is explicitly taken into 
account, with the result that a constant fallin g/ascending velocity sets in, th us dom- 
inating the particle long-time transport (iMaxev fc CorrsinI 1986: . Maxev 1987alftl. 1990 ; 



20071: iMartins Afonsoll2008l) . 



Aliseda et atll2002l:lFriedman fc Katdl2002l : lRuiz Macias fc Peter£i2004l:lMarchioH Fantoni fc Soldati 



Our main aim here is to focus on transport of inertial particles in compressible media. 
Although our results are rather general (i.e. they apply to both one, two and three di- 
mensions, and for a variety of external forces), the present study is physically motivated 
by the transport of floaters on a surface flow (e.g. the ocean surface) in the presence of 
strong surface winds. As we will see, the resulting dynamical equations for the particles 
moving on the (horizontal) surface under the action of strong (constant) winds are for- 
mally identical to those of heavy particles moving along the vertical under the action of 
gravity. A constant drift is thus expected, analogous to the settling velocity of a particle 
in suspension, in the presence of gravity. 

The specific problem we aim at investigating can be summarized as follows. The drift 
velocity of a floater on a still water surface will be in general a complicated function 
of the wind velocity, (7, and of other parameters like e.g. the floater structure and the 
surface roughness. A similar complication is expected in the way in which the drift ve- 
locity adapts to the wind and surface current variations. A rather minimal model could 
be obtained assuming a simple relaxation dynamics for the drift velocity: 

v^[V-{v-u)]-f, (1.1) 

where u is the water surface velocity, 7 = 7(C/ — u,v — u) is the relaxation rate, and V 
is the terminal drift velocity (with V — V{U — u)). Since we expect V,u<^ U, we may 
approximate J7 — it ~ [/ in the arguments of both V and 7, and, as a rough approxima- 
tion, we may consider a linear relaxation dynamics: j{U, v — u) ~ "y{U, 0). Under these 
hypotheses, p.ip becomes formally identical to the one for a small heavy particle in a 
viscous fluid, in the presence of a gr avitational acceleration 7V. In this case, 'y~^ would 



be the Stokes time (jMaxev fc Rilevi 1983: Michaehdes 1997). 



Now, spatio-temporal variations in the surface water velocity will lead to difficulties in the 



determination of a mean drift velocity, due to preferential concentration effects (jMaxev 



1987tl) . In other words, the mean drift may not necessarily be equal to what would 



be obtained by calculating a spatio-temporal average. In general, the interplay between 
currents and particle trajectories might lead to the result that either an enhanced or a 
reduced drift (with respect to the one in still fluid) might appear. 

Here, we will assess the abo ve possi b ility by means of an analytical (perturbative) ap- 
proach in the same spirit of iMaxevI ( 1987a ). We will be able to identify an important 



dynamical feature of the carrying flow (the way through which it decorrelates in time) 
responsible of the different behaviour of the resulting drift with respect to the corre- 
sponding value in still fluid. 

The paper is organized as follows. In §[2] the basic equations governing the time evo- 
lution of inertial particles in a prescribed flow are given, together with the perturbative 
expansion for strong sweeping or gravity. In § |3] we focus on compressible flows and 
compute the leading correction to terminal velocity. The same quantity is calculated in 
§ |4] for incompressible Gaussian flows. In § [5] we analyse the phenomenon of particle 
diffusivity and, at the leading order, we study the diffusion coefficient; we also perform 
a quantitative comparison of the drift and diffusion displacements. Conclusions follow in 
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§ [6l The appendices are devoted (§ |A| to displaying some additional analytical results, 
and (§ El, to provide some technical details of the calculation. 



2. General equations 

In the hypotheses of (|l.ip . the motion of a floater dragged by the wind, on a water 
surface with velocity field u{x,t), will be described by: 

x{t)=v{t) , 
v{t)=l[V + u{x{t),t)^v{t)], ^'■'> 

with the quantities V and 7 assumed as constants. We shall denote the direction of the 
vector V as "sweeping" , and we align the axes such that it corresponds to the positive 
Xd component. To make contact with the dynamics of a heavy particle in a viscous fluid, 
we maintain d, that is the number of dimensions, arbitrary. 

Given (Eulerian) characteristic length and time scales L and T and amplitude cr„ for the 
velocity field tt, we introduce the dimensionless Stokes, Kubo and Froude numbers 5, K 
and F, defined as 

5=—, and F = . (2.2) 

7T' L VtV^ 

Notice that such a definition of F is consistent with the "identification" g ~ 7V, if 
gravity is the "sweeping" force under consideration (rather than wind). 

From now on, we adimensionalize times with 7"^ and velocities with (7„, and denote 
the new adimensional variables with the same letters as before. In dimensionless form, 
the characteristic scales L and T of the velocity field, and the bare terminal velocity V, 
will read: 

L = S-^K-'^, T^S-^ and V ^ SKF'^ . (2.3) 

The Kubo number K is basically the ratio of the life time and rotation time of a vortex, 
with K Q corresponding to an uncorrelated regime, K —t 00 to 'a frozen-like regime, 
and real turbulence being realized by K — 0(1). 

We confine ourselves to situations in which the particle trajectory, after some transient 
regime depending on the initial conditions, shows small deviations from the "sweeping" 
line, thus allow ing us to deal with a quasi- ID problem as a zeroth-order approximation 



( Maxev 1987a ). This happens when the effect of streaming (or gravity) is much stronger 



than the deflections due to the external flow. 

We thus isolate in the solution v = v{t) of (|2.ip a term associated with the deviation 
from the behaviour in still fluid: 

i(t) = f dt'e''''u{x^°\t') + xit'),t') , (2.4) 

where 

a;""(i) = a;(0) +Vt+ [v[0) - V](l - e"*) (2.5) 
accounts for the unrenormalized sweep, and Ax/dt = {j, so that: 

x{t)^[ dt' il^{t-t')u{x^°\t') + x{t'),t'), VW = l-e"*. (2.6) 

A perturbative solution of (|2.6p rests on the smallness of x{t) in the argument of u. 
More precisely, it is necessary that i(rp) <C i, with Tp the correlation time for the fluid 
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velocity u{x{t),t) sampled by the particles. We can estimate 

Tp ^ min(r, L/au,%^), %^ = L/V = S-^R-^F^ 



(2.7) 



with giving the contribution from sweep to decorrelation. In the absence of sweep 
effects, we know that x(t-o]/L will be small provided eit her if <C 1, or JsT > 1 with 
SR-'^ > 1 (jWilkinson fc Mehlidliool lOlla fc Vuoioll2007[). I n both regimes, indeed, the 
inertial particles will see m as a Kraichnan field ([Kraichnan Il968lll994l) . Sweeping acts 
to reduce the correlation time. A sufficient condition for small x(Tp)/L turns out to be, 
in this case: 

SRF~^ > 1 , (2.8) 

that is the strong-sweep condition V ^ (t„ (no condition is put on 7). Notice that, for 
sufficiently large S, this condition and the one for a Kraichnan regime JsT <C 1 overlap. 
We assume m as a homogeneous, isotropic, stationary, zero-mean random flow and we 
denote with (•) the ensemble average over its realizations. We are interested in finding 
the steady-state average particle velocity, which corresponds, in (|2.4p . to averaging over 
u and taking t +00. 

Taylor expanding in x the right-hand side (RHS) of (|2.4p . and using (|2.6p recursively, 
allows us to calculate the correction to sweep. We clearly obtain (^'^'(t)) — 0. Then, we 
have 

(«rW) = f dt'i^{t-t'){U,{t')d,Um) , (2.9) 

where 

U{t) ^u{x^°\t),t) (2.10) 

(i.e. the unperturbed flow is "sampled" on the fixed "sweeping" line). As it is well known, 
the lowest order correction to V vanishes if u is incompressible. In this case, in order to 
obtain non-zero corrections to the falling v elocity, it is ne cessary to go to higher orders in 
the perturbative expansion of (ISlll-jMl) (|Maxevlll987^ . Namely, in the incompressible 
case, we have: 



{vf (t)) = / At'i,{t-t') dt" m' ^t")-i,{t~ t")] {U'^U'MU,) 



(2.11) 



The latter quantity is zero if m is a Gaussian field. Thus, for incompressible Gaussian 
flows, one has to compute the next order: 



= /*dt' f dt" i,{t~t')ij{t' -t") X 

dt"' i^it - t"') - f dt"'^{t" - t"') 
10 Jo 



(2.12) 



X 



In §[3] we provide an example of compressible flow, for which (|2.9p applies. Due to the 
difficulty of dealing analytically with non-Gaussian flows, we will not provide applications 
of (|2.1ip . On the contrary, in § we will focus on an incompressible Gaussian flow, for 
which (|2.12p is the leading correction to the falling velocity. 



3. Compressible flows 

Let us focus on the compressible case and compute the average in (|2.9p . Clearly we 
only need to compute the component with index i = d (in our convention it is assumed 
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positive if pointing along the mean flow). 

Let us first analyse d > 1. Introducing the well-known compressibility degree, V = 
{{djUj)"^) / {{dkUi){dkUi)) e [0,1], the expression for the Eulerian pair correlation tensor 
Rij{x,t) = {ui{x,t)uj (0,0)) can be deduced from 

R^J{x,t) = [(1 - dV)^^^J - (1 - V)S.,jd^] n{x,t) (3.1) 

by imposing the form of the scalar TZ{x, t). Let us assume a Gaussian behaviour both in 
space and in time, with temporal oscillations described by uj (adimensionalized with 7 
for the sake of consistency) : 

Substituting into (|3.ip . we obtain 

i?.,(^,t) = ^^^cosMe-^^*V2e-^^^^-^/2 X (3.3) 
X {(1 - dV)S^K^XiXj + % [{d - 1) - (1 - V)S^K^x^] ] 

and thus 

d,R,,{x,t) = ^^!^cos(^i)e-«'*'/2e-^'^'-'/2 [^2 ^2^2 _ (^^3)] x, . (3.4) 

The case d — 1 automatically implies V = 1, and expressions p.ip through (|3.3p are ill 
posed. However, if one neglects them and considers expression (|3.4p directly, everything 
is consistent and also the ID case can be investigated by means of the same formalism. 
Our choice of obtaining the latter equation passing through the former three for d > 1 
is simply dictated by the simple and important meaning that the compressibility degree 
plays: the lower bound V = Q denotes incompressible flows and the upper bound V = I 
perfectly compressible (potential) ones. 

It is interesting to study under what conditions the tensor p.3p is positive definite (for 
d > 1). One can notice that, as in our calculations it always appears under some integral, 
it is sufficient to study the separations x < L, because the Gaussian factor makes the 
contribution from larger distances neghgible. Thus, using (|2.3p . S'^K^x^\x<l < 1, and 
therefore the quantity in square braces beside Sij is always positive in this case. Then, 
the investigation of the tensorial part suggests that Rij would surely be positive definite 
(for the only relevant interval, x < L) when V < d~^, if the time-oscillating factor 
were absent. This tells us that, at low compressibility, the introduction of the temporal 
oscillations corresponds to the appearance of negatively-correlated areas, and tuning the 
parameter lu one can change the relevance of the recirculating flow. 

Under our assumptions, performing the simple change of variables t = t — t' , at very 
large times Vd expression (|2.9p for the "sweeping" component simplifies to 



where 



F4 



- (d + 2) 



cos(wT)V'(T)e-^'"'/2 , (3.5) 



F = S^Jl + A2 w and A^SK^p-"^ (3.6) 

are respectively the decorrelation rate, and the ratio T/%w 

The last integral in (|3.5p can be carried out analytically and the result can be expressed 
in terms of error functions, but for the sake of simplicity it is not reported here (see (jA 1| 
in appendix |A| . Out of the six parameters in play, we notice that V is exactly linear in 



6 M. Martins Afonso, A. Mazzino and P. Olla 
0.003 I 1 1 1 




-0.005 



2 4 6 8 10 

Figure 1. Leading-order correction to the terminal velocity, normalized with the 
corresponding bare value, as a function of u, for d = 1, 2 and 3. Here, V = 1 = S and 

F = 0.1 = K. 



V , and coherently vanishes in the incompressible case, thus we can fix P = 1 for the rest 
of this section without loss of generality. Let us then analyse the dependence of V on the 
five remaining parameters {uj, S, K, F and d). 

• A simple look at the structure of the integrand in (|3.5p . or at the final expression 
(|A ip . suggests that V ^ tor uj oo. Equation (jA2| displays the simplification 
occurring in the static case a; = 0. In figures [T] and [21 we plot the ratio V/V between 
the leading correction to the terminal velocity and its bare value, as a function of the 
frequency. We notice that small values of the frequency are associated with a negative 
renormalization, while high-frequency flows are characterized by an increased terminal 
velocity. A similar situation can be shown to arise as a result of finite-T correction s, in 
the transport of inertial particles by a Kraichnan velocity field (jOlla fc Vuoldl2007l ). We 
can thus conclude that (at least for P < d~^), if the extension of negative regions in the 
time correlation profile is sufficiently small (large), the renormalized streaming is smaller 
(larger) than the bare one, which amounts to saying that the absence (presence) of areas 
of recirculation helps to make the particles travel slower (faster) than in the case of still 
fluid. As we can see from figure [H this conclusion is confirmed also for finite values of 
the expansion parameter; in the case in exam: SKF~^ = 4. Notice that an oscillating 
exponential profile for the time correlation has been utilized in this case, in place of the 
oscillating Gaussian of (|3.2p . 

• For T ^ TsMj (see (|2.7p ) the integrand in (13. 5p will depend solely on T^w Thus, if 
one studies the zero of expression (i.e. the boundary that separates accelerating and 
decelerating situations, at the lowest order), the most interesting plot is represented by 
the separatrix line in the plane lo vs. SKF~^ for not-too-small values of K, such that 5*, 
K and F appear only via this combination. In figured the lower right area corresponds to 
deceleration and the upper left to acceleration in sweeping. Notice that low (respectively, 
high) values of SKF~^ correspond preferentially to situations of increased (respectively, 
decreased) terminal velocity. 
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Figure 2. Same as in figure [T] but with S = W and F = 1 = K . 
0.25 



0.2 



0.1 
^ 0.05 


-0.05 
-0.15 



1 


1 1 


1 


1 1 

^ < 


















- ^ 
> 

1 1 1 1 1 1 1 



0.5 



1.5 



2 

a; 



2.5 



3.5 



Figure 3. Leading-order correction to the terminal velocity vs. ui, from a numerical simulation 
of transport by a 2D fully compressible velocity field. Parameters of the simulation: 5 = 2, 
K = 0.5 — F. The size of the spatial domain is 10 L. 



• The same simplification into a dependence on SKF~^ only is not valid for the whole 
expression p.5|) itself, because of the prefactor. Nevertheless, an interesting rescaling 
property can be shown to hold when the parameters belong to a specific range. First of 
all, if one considers the ratio V/V, the prefactor S^K^F~^ in p.Sp becomes S^K^ , i.e. 
{SKF~^)^F^. Now, suppose that T ^ 7^^, such that, in the exponent of the Gaussian 
factor, A > 1 and thus F ~ {SKF"^)"^. Also, suppose that SKF"^ > 1, so that 
the same Gaussian is very narrow. This allows us to Taylor expand ^{t) ~ r. Then, 
starting from a situation in which such approximations hold (e.g., the case plotted in 
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Figure 4. Separatrix between regions of increased (upper left) and decreased (lower right) 
sweeping in the plane to vs. SKF~^ , for d = 1, 2 and 3 (solid, dashed and dotted lines, respec- 
tively). The plot namely corresponds to fixed values of K — 1 and F = 0.1, but varying these 
parameters only negligible changes appear (unless K <^ 1). 

figure [2|), suppose to rescale SKF~^ by a factor a. The expression in square brackets 
and the simplified Gaussian factor remain unchanged upon rescaling t by a~^, and 
consequently the cosine does not change if uj is rescaled by . Keeping into account the 
T factor exphcitly appearing in (j3.5p . the one coming from the Taylor expansion of ipir), 
and the integration variable differential, one deduces that the integral gets rescaled by 
{a~^)^ = a~^. By considering also the above-mentioned prefactor {SKF~^)'^F'^ , one can 
argue that, under these approximations: I) upon rescaling SKF~^ by a while keeping 
F fixed, figure m simply reproduces itself by rescaling u by on the abscissae and V/V 
by a"** on the ordinates; II) upon rescaling F by P while keeping SKF^^ fixed, figure [5] 
simply expands vertically by without changing horizontally. 

In figures [5] and m V/V is now plotted as a function of SKF~^ for different values of the 
frequency. 

• The dependence on the dimension seems to be very weak. In particular, all our plots 
show little difference between d = 1, 2 and 3, with the one-dimensional case showing 
the largest peaks. Moreover, figures [THU and [5HB] are peculiar in the sense that the three 
lines always cross in the very same points. We have not been able yet to identify the 
underlying analytical mechanism causing this phenomenon, which however seems to be 
very robust because it is also confirmed by plotting e.g. a fictitious four-dimensional case. 



4. Incompressible Gaussian flows 

The role of w, i.e. of oscillations in the velocity correlation profile, in determining 
the sign of the terminal velocity renormalizatio n, is completely different in the case of 
a Gaussian incompressible flow. As discussed in iMaxevI (|l987n . this requires going to 
fourth order in the perturbation expansion (see (|2.12p ). 

The interest of the problem is rather formal, as non-Gaussianity in realistic turbulent 
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Figure 5. Leading-order correction to the terminal velocity, normalized with the corresponding 
bare value, as a function of SKF~^ , for d — 1,2 and 3. Here, P = 1 and uj = 10. The plot namely 
corresponds to fixed values of = 0.1 — F, and varying these parameters it gets rescaled as 
explained in the text. 
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Figure 6. Same as in figure [S]but for the static case a; = 0. 



flows is likely to play an important role. Nevertheless, simple analytical expressions can be 
obtained in the small-sweep regime, with rapid decorrelation imposed through smallness 
of T. More precisely, we assume 



SK^P-^ < 1, 5 > max(l, K^) 



(4.1) 
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Figure 7. The profile of the function Q. 



In this limit, the decorrelation T is dominated by the Eulerian correlation time T <C Tg = 
1 and the total displacement in a time T will be much smaller than the characteristic 
scale L = S-^R-^. 

As detailed in appendix|Bl this allows to expand the integrand in (I2.12p in both space (the 
correlations {UU)) and time (the propagator -0). Notice that, in these approximations, the 
precise form of the spatial correlation (jB ip has no influence on the sign of the correction, 
as only its second derivative computed at merged points, (|B 3|) . enters the calculation. 
A simpler expression for the setthng velocity renormahzation, illustrating the changing 
sign behaviour, is obtained with a temporal correlation in the form: 



hit) 



- ^-s\t\ 



cos(wi) 



Substituting into (|B 7p . we obtain 



V= lim 



\t)) 



{d+l){d + 2)f^K^ 
4(d- 1)S'2F2 



Q 



where 



Q(a) 



11 - 94a2 + 177a4 - 450^ 



7a« 



UJ 

S 

,10 



(4.2) 



(4.3) 



(l + a2)7 

Using (|2.3p . we see that V /V cx S^'^K'^, corresponding to an 0{S^^ K*) renormahzation 
of 7, to be co mpared with the 0( S~^K^) renormahzation observed in the fully compress- 
ible ID case ( Olla fc Vuolol 2007 ). The plot of the function Q{a) is shown in figure [71 It 
was observed in iMaxevI (|1987^ in the case of strong sweep, and in the absence of oscil- 
lations for h{t), that the correction to the settling velocity is positive. We see in figure 
[7] that the same situation occurs in the small sweep regime, while we can also identify a 
range of values of uj for which the falling velocity is decreased. However, the renormahza- 
tion effect quickly becomes negligible at large cu. In any case, we were able to show that, 
at least in some particular instances, the presence or absence of compressibility makes 
the role of ui the opposite passing from one case to the other. 
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5. Particle diffusivity 

After computing the terminal velocity, it is interesting to study how the particle diffuse. 
In particular, we would like to find the diffusion coefficient 



D 



lim 

t — >co 



2dt 



(5.1) 



to be compared with the diffusivity of a fluid parcel, 

V ^ alT = S-\ 

With this aim, we firstly notice that the term a;<°> in (|2.5[) does not contribute to expres- 
sion (|5.1[) . because it is independent of the external random flow. It is thus sufficient to 
investigate x from (|2.6p . and namely (using (|2.10p ') its leading order 

x^^>{t) = f dt' iP{t - t')U{t') . 
Jo 

For zero-mean flows we have (i<^'(i)) = 0, therefore: 

i^d) ^ lim -L f At'4'{t-t') f dt"i^{t-t"){U,{t')U,{t")) . (5.2) 

Looking at (|5.2p . we see that ijjit — t') and ip{t — t") differ from 1 only for t' and t" close to 
and for t — > oo we can disregard this contribution. Thus, D'"^' is simply the diffusivity 
of a random walker with velocity lY, and we have the standard expression 



1 



dt{umu^{t)) 



(5.3) 



If we now restore our space-and-time Gaussian oscillating correlation, we can replace the 
average in (|5.3p with the trace of expression (|3.3p . computed at space-time separation 
iyt.t), and we obtain the result: 



V2d2r5 



(2d - 1) 



UJ\2 

S 



(5.4) 



with r and A given in p.6p . 

Notice that the compressibility degree V disappears in the trace Ru, and thus does not 
affect the diffusion coefficient. It can easily be seen from ()5.4|) that I?'^' depends on S and 
SKF~^ only, and not on K and F separately. Figures [8] and [9] show the behaviour of the 
diffusivity as a function of the frequency. In figures [10] and [Til the diffusivity is plotted 
versus SKF^^. In all cases, plotted is the ratio of the particle-to-fluid diffusivities. 

As a further step, one can also investigate the different diffusive behaviour between 
the "sweeping" direction and the orthogonal one(s) (if any). In particular, one obtains 



= lim 



it)] 



and 



2dt 



V2d^T5 



{2-V)+V{- 



1 



V2d^T5 



/2r^ 



{d + V-2)A'' 



{2d + V-3) + {l-V){^y 



A2 + (d-l) 
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Figure 8. Leading-order expression for the total diffusivity, divided by the corresponding 
fluid-particle value, as a function of u, for d = 1, 2 and 3. Here, S = 10, — F and VP. 




(which correctly vanishes for d = 1 because there V — 1). 

It is worth noticing that the compressibihty degree now appears again in the two previous 
expressions: only their sum is independent of it. 
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Figure 10. Leading-order expression for the total diffusivity, divided by the corresponding 
fluid-particle value, as a function of SKF'^, for d = 1, 2 and 3. Here, = 10 and 5 = 1, VP. 




Figure 11. Same as in figure [TU] but for the static case a; = 0. 



5.1. Comparison between drift and diffusion displacements 

It is now possible to compare the contributions to a particle displacement given by drift 
and by diffusion. At the time t large enough to forget about initial conditions, the typical 
displacement due to drift is about (V -I- V)t, while the diffusive one is of the order of 
\/2dD'-^H. The former clearly dominates at large times, while the latter represents the 



14 M. Martins Afonso, A. Mazzino and P. Olla 




2 4 6 8 10 

u) 

Figure 12. Crossover time, between diffusion- and drift-dominated particle displacement, as a 
function of ui. Here, d = 2, V = I, F = Q.l and SKF'^ = 1. 

leading value of the displacement at short times. We can thus define the crossover time 

" (V + vy ' 

after which drift represents the dominant mechanism and diffusion can be neglected. 
Focusing on the two-dimensional potential-flow case, in figure [T^] such crossover time 
(remember that, in our convention, it is adimensionalized with 7) is plotted as a function 
of the frequency for three different values of the Stokes number. Notice the peak that 
can occur in t, at a critical value of ld when the other parameters assume certain values. 



6. Conclusions and perspectives 

We investigated the drift and diffusion behaviour of inertial particles in random flows, 
in the limit in which the expansion "a la Maxey" works, i.e. in the presence of strong 
sweeping, or of an underlying flow varying rapidly, such that the resulting motion at the 
scales of interest is almost rectilinear. This picture can be applied both to the settling of 
heavy particles due to gravity, and, more interestingly, to the motion of floaters, e.g. on 
a water surface, under the action of a streaming wind. 

Imposing an oscillating Gaussian form to the external flow, we found the analytical 
expressions of the renormalized terminal velocity and of the effective diffusivity, as func- 
tions of the oscillation frequency and of the Stokes, Kubo and Froude numbers. The 
latter three quantities often appear not independently from one another, but rather as 
an overall combination. For the terminal velocity, we found that, if the flow is compress- 
ible, small values of the frequency are associated with a negative renormalization with 
respect to the corresponding value in still fluids, while high-frequency flows are character- 
ized by an increased asymptotic speed. In some cases, the opposite behaviour holds as a 
function of the parameter SKF^^, or as a function of ui itself if the flow is Gaussian and 



Renormalized transport of inertial particles in surface flows 



15 



incompressible. The effective diffusivity can show some "resonant" peaks as a function 
both of oj and of SKF~^. The role played by the compressibility degree V is apparently 
quite simple: it appears linearly in the leading correction to the renormalized terminal 
velocity (so that, for incompressible flows, higher-order effects have to considered), and 
it does not influence the total effective diffusivity, because the dependence on V shown 
by the components of the diffusivity parallel and orthogonal to the sweeping direction 
is such as to create an overall compensation. However, we were able to show that the 
value of V can have relevant consequences on the role played by to on the renormalized 
streaming. In particular, we can conclude that, for compressible flows (at least when the 
compressibility degree is not too high, but neither too low in order to have non-negligible 
effects), the introduction of time oscillations in the velocity correlation function mimics 
the presence of areas of recirculation, or in other words of negatively-correlated regions 
in the flow, and is associated with an increase in the particle terminal velocity; on the 
contrary, a decrease is found for slow or no oscillations (for which the Eulerian time 
correlation is positive). In all cases, we also found that the dependence on the dimension 
{d = 1, 2, 3) is quite weak. 

MMA is grateful to Keck Foundation for financial support and to Grisha Falkovich for 
useful discussions. 



Appendix A. Exact results for the compressible case 

The resuh of the integral in (|^ . computed with MATHEMATICA (which makes use 
of complex variables), is: 



dF^r^ 



V 



{d - 1) + ({2d +1)+ (i^) A'' + {d + 2) 



x(l-iw)e(i-'-)'/2r^ h - erf 



1 — icij 
V2T 



-FA' 



^^2^e-"'/2r' X 



(d-l)A4+ [ (2d+l)+ (|)MA2 + (d + 2) 



erf 



(v2r) 



(Al) 



where erf is the usual error function, 3? denotes the real part, and F, A were defined in 
(ESI). 

In the static case (t^; = 0), this expression simplifies to: 



1 — erf 



1 



V2T 



(A 2) 



(rf - 1)A^ + ( (2rf + 1) - — ) A^ + (d + 2) 



Appendix B. Technical details for the incompressible case 

Setting V ^ 0, we can rephrase (|3.ip following iMaxevI ( 1987^ and write 



Ri 



^ix,t) = [g{x)6,, + [fix) - g(x)]^}hit) , 



(Bl) 
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where f{x) and g{x) are such as to ensure incompressibihty, and we assume h[Q) = 1. In 
this way, the integrand of (I2.12|) . with i — d, can be rewritten as: 

(B2) 



H{t"',t\t",t) = dMW)d'k{u'jUn 

= sgn(<' - t"')h{t - t")h{t' - t'") 



r a — 1 



with s — V\t' — t"'\ and r — V\t — t"\ (primes on / and g now indicate differentiation 
with respect to the corresponding variables). Notice that H is invariant under changes of 
third and fourth variables, while it flips its sign when changing first and second variables; 
homogeneity also implies that H{t"', t' , t", t) is a function of t — t" and t' — t'" only. 
The first condition assumed in (|4.ip . corresponding to a quasi-uniform field, allows us to 
expand the functions / and g in (jB 1\ in powers of the spatial variable: 



+ j[/"(o) - 9"mnT, 



(here, /'(O) = = g'{0) to ensure smoothness, and /(O) — g{0) — d^^ by construction). 
The incompressibihty constraint gives us g"{0) — [{d + l)/{d— l)]/"(0). Reminding the 
meaning of L and our adimensionalization of lengths in (|2.3p . we can introduce 

/ = lV"(0) ^ f"{Q)^S^K^f, (B3) 

where / is now 0(1). 

Looking at (jB 2[1 . we see that we must expand the first-order derivatives; this gives us: 



H{t"', t\ t", t) = {t' - t"')Vh{t - t")h{t' - t'") 
{d+l){d + 2)S^K'^P 



(d + i)/"(0)g"(0) 



{d-l)F^ 



h{t-t")h{t' -t"'){t' -t'") 



(B4) 



We then see that the contribution from the space correlation has always the same sign. 
The integrals in (|2.12p can be written in the equivalent form: 



{if\t))^l drV'(r) 
At" 



At' 



At" iI}(t + t')iI){t" -t-t') 



At' V^t" - t)^{t + t' - t") \ H{t', 0, 0, r") 



leading to the expression for V = \imt^oo{vd (t)): 



V 



At V'(t 
At" 



At' 



At" ^(r + T')i:{T" ^ t - t') 



At' 7/;(r" - r)7/;(r + r' - r") \ H{t', 0, 0, r") 



We can eliminate the integral in r rewriting (see figure [T3|) : 



At / At' 



At" = 



r-\-T' 



At" 



At' 



At 



max(0, — r') 



and 



At / At" 



At' = 



At" 



At' 



At 



max(0,T''' — r') 



(B5) 



(B6) 
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Figure 13. The domain for ([B 5|l is the wedge between planes A and B and the coordinate plane 
t't" . The domain for (|B 6|l is the wedge between planes B and C and again the coordinate plane 
t't" . The intersection of both planes B and C with the coordinate plane tt" is the line t" = t. 



Carrying on the integrals in r by exploiting (|B 4[) : 

^ (l+lKl+^l^!:^ Td ''J 'r* r ' 

V = ^^^I)]?^ X Mr )|y^drG,(r,r) 

+ dT'[G2{T',T")+G4{T\T")]+J dT'G3{T',T")jT'h(-T') , (B 7) 

where the kernels (/i = 1 , . . . , 4) get the following, simple expressions as the second 
condition in (|4.ip allows us to Taylor expand the propagators ip in r' and r": 



Gi(r',r") 
G2(r',r") 
G3(r',r") 
G4(r',r") 



— r' 
r '— r 



drV'(r)-0(r + t')%1>{t" - t -t') 
dTi){T)i){T + T')V'(r" - r - r') 
dr V'(t)V'(t" - r)i/>(T - t" + r') ~ 







6 






t"V" 






T2 + 


6 



'3 ^,/4 



6 12 



r T r 



6 12 

r" 

dr V-MV^r" - t)?/'(t - r" + t') ~ — + 



// 

r T T 



12 6 



The desired result can then be obtained from (|B 7p by specifying the exact form of the 
temporal correlation /i; imposing (|4.2p . we finally get (|4.3p . 
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